NPS ARCHIVE 
1959 

ST. VILLE, E. 






h.l Na I 

M ■ 












I 






854 



HEAT TRANSFER CHARACTERISTICS 
OF SCREEN MATRICES AT LOT REYNOLDS NUMBERS 



*A” 



Edward L. St. Ville 



f 



r 




h i 1 



ll 



* ' Y L • 




HSAT TRANSFER CHARACTERISTICS 
OF SCREEN MATRICES AT LOV, T REYNOLDS NUMBERS 



by 



Edward L. St. Ville 

n 

Lieutenant, United States Navy 



Submitted in partial fulfillment of 
the requirements for the degree of 



RASTER OF SCIENCE 
IN 

MECHANICAL ENGINEERING 



United States Naval Postgraduate School 
Monterey, California 



19 5 9 



HEAT TRANSFER CHARACTERISTICS 
OF SCREEN MATRICES AT LOW REYNOLDS NUMBERS 

by 

Edward L. St. Ville 

This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 

MECHANICAL ENGINEERING 
from the 

United States Naval Postgraduate School 



ABSTRACT 



The convective heat transfer characteristics of six screen matrices 
with a range of porosity from 0.602 to 0.832 are presented. The charac- 
teristics were investigated over a Reynolds number range of 0.1 to 10, 
based on hydraulic diameter. The data was obtained by passing air 
through a 6 in. diameter woven screen matrix composed of three stacked 
screens. 

The heat transfer data was obtained by use of the transient techni- 
que developed by G. L. Locke at Stanford University. Temperature-time 
history of the downstream air temperature is used to determine the maxi- 

i 

mum slope of the experimental cooling curve. From a plot of maximum 
slope versus NTU of the generalized cooling curve, which is a unique 
solution for this type of cooling, an evaluation of the experimental NTU 
is possible. With the NTU thus determined, the experimental value of 
the thermal conductance coefficient, h, is readily obtained. 

The present data is an extension of existing data and was obtained 
using the same type of matrices but with a slightly different technique 
and over a lower range of Reynolds numbers. The present technique in- 
volves heating the screen matrix to a uniform temperature and introduc- 
ing the screens into a stream of air entering at a uniformly lower tem- 
perature. The temperature-time history of the downstream air is recorded 
and available for analysis. 

Application of the present data is thought to be useful in transpira- 
tion cooling of environmental chambers subjected to an external heat 
source. 
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TABLE GF SYMBOLS AND ABBREVIATIONS 



English Letter Symbols 

A - heat transfer area, ft* 

Af matrix free flow area, Af = a. A^, ft^ 

Aj’ - free flow area of a single screen or a perfectly packed 
screen matrix, Ajp’=. cc A^, ft 2 

A^ - matrix frontal area, ft 2 

Cf - specific heat at constant pressure, BUJ/(lbm °F) 

c 8 - specific heat of solid phase, BTU/lbm °F 

C - thermal capacity rate, C = wCp, BTU/(hr °F) 

D - Diameter, ft 

D^ - hydraulic diameter, D^ = A r^, ft 

Dh* - hydraulic diameter of single screen or of a perfectly 
packed screen matrix, D^’ = 4 r^’ , ft 

G - mass velocity, lhm/(hr ft^) 

G^ - mass velocity based on matrix free flow area, 

Gi - Wf/Af, lbra/(hr ft 2 ) 

G^’ - mass velocity based on free flow area of a single screen 

or of a perfectly packed screen matrix, 

Gi’ = W f/Af * , lbra/(hr ft 2 ) 

h - unit conductance for convective heat transfer, BHl/(hr ft 2 °F) 
k - thermal conductivity, BTU/(hr ft 2 °F/ft) 

L - flow length of matrix, ft 

n - number of screens in matrix 

q - heat transfer rate, BTU/hr 

r h “ hydraulic radius, r^ = or//? , ft 

r h* “ hydraulic radius of a single screen or of a perfectly 
packed screen matrix, r h ' = ^ '/ ^ * , ft 

t - temperature, °F 

T - temperature, °R 
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V - velocity of fluid based on frontal area, ft/sec 

- velocity of fluid based on free flow area, ft/sec 

Vj_* - velocity of fluid based on free flow area of a single 

screen or of a perfectly packed screen matrix, ft/sec 

vtf - mass flow rate of fluid, lbra/hr 

VI f - mass of fluid contained in pores of matrix, lfcro 

W s - mass of solid in matrix, lbm 

x - distance from entrance of matrix to a given cross section 
measured in the direction of flow, ft 

Greek Letter Symbols 

C( - porosity, oc = volume of voids/ total matrix volume, nondlmen- 
sional 

cc‘ - porosity of a single screen or of a perfectly packed screen 
matrix, nondime ns ional 

j3 - heat transfer area per unit volume of matrix, ft”-*- 

ft' - heat transfer area per unit volume of a single screen of or a 

perfectly packed screen matrix, ft”-*- 

Q - time, hrs 

- viscosity, lbm/ (hr ft) 

- mass density, lbm/ ft^ 

Dimensionless Groupings 

NTU - number of heat transfer units, NTU = hA/Cw^c^) 

Np r _ Prandtl number, Npj. = Cf y^f Af 
%e - Reynolds number, % e = 

Uro* - Reynolds number employed for screen matrices, 

Ng t - Stanton number, N s ^ = h/GjCf 

%t* “ Stanton number employed for screen matrices, Ngt' = h/Gj/Cf 
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T - generalized time variable, t' = hA/W s c 3 ( 9 - Wfx/w-fL) 
z - generalized position variable, z = (hA/w f c<0 x/L 
Subscripts 

a - refers to atmospheric conditions 

f - refers to the fluid phase 

i - refers to matrix initial state 

s - refers to the solid phase 

1 - refers to the fluid state upstream of the matrix 

2 - refers to the fluid state downstream of the matrix 



i 



1, Introduction 



In this study, porous media is defined as a solid matrix possessing 

2 3 

a high ratio of void surface area to bulk volume ft / ft , and a void 
geometry which permits flow of gas or liquid through the matrix. The 
term porous media is used throughout, to refer only to those random homo- 
genous configurations which result in three dimensional flows, with fre- 
quent contractions and expansions. 

In the past, heat transfer characteristics of porous media have been 
of interest primarily in the study of the flow of oil and water through 
underground formations, the flow of gases through blast furnace charges, 
and the flow through the granular packing used in absorption and stripping 
columns and catalyst beds. Other applications of more recent interest are 
found in; 

a. Periodic-flow, regenerative type heat exchangers 

b. Nuclear reactors 

c. Boundary layer control for air foils 

d. Transpiration or sweat cooling of environmental enclosures. 

In regenerative heat exchangers the energy transfer from the hot to 

cold fluid is accomplished first by transfer from the hot fluid to the 
solid matrix followed later by a transfer from the matrix to a cold fluid. 
Thus, the matrix is alternately heated and cooled either by periodic move- 
ment of the solid from the hot to cold streams and back again, as in 
rotary- type regenerative exchangers (Ljunstrom air-preheaters), or by 
sequential valving of the hot and cold gas streams through two solid 
matrices in parallel. In the past, types of matrices used were of the 
nature of brick checkerwork or packed corrugated plates — constructions 
which do not possess a very large transfer area per unit volume. As an 
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illustration of the wide variance in the degree of compactness of various 
matrices the following tabulation was prepared. 



APPROXIMATE 
AREA/UNIT VOLUME 



MATRIX ft 2 /ft 3 

1. Brick Checkerwork 0.8 

2. Rod banks, rod dia. - 0.375 in. 51 

* = 1 = 2 

3. Corrugated plates 150 

(Used in Ljungstrom heaters) 

4. Sphere Beds 

a. 1/4 in. spheres 170 

b. 1/16 in. spheres 690 

c. 1/32 in. spheres 1380 

5c Packed Wire Screens 

a. 5 mesh 190 

b. 20 mesh 760 

c. 50 mesh 1900 

d. 300 mesh 11,400 



- Ratio of longitudinal pitch to tube diameter for flow 
across banks of circular tubes, dimensionless. 

- Ratio of transverse pitch to tube diameter for flow across 
banks of circular tubes, dimensionless. 

Nuclear reactors for electric power generation consist essentially 
of a solid matrix in which the thermal energy created by fission is trans- 
ferred as heat to a working substance passing through it. Thus the possi- 
bility of making the fuel in the form of a porous medium through which to 

pass the cooling fluid should be considered. The principle advantage 

0 

being that of high heat transfer effectiveness which would enable the use 
of a gaseous medium for a primary coolant. In addition, if porous solids 
are practical for this application it may be possible to reduce the volume 
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and the overall bulk and masa of the system. 

The feasibility of using porous media for boundary layer control 
has come to the fare primarily due to the development of sintered metals. 
Sintered metals are an excellent example of porous media, the area density 
being of the order of the 300 mesh screens. Interest has centered princi- 
pally in the use of area suction to delay boundary layer transition and 
thus take advantage of the lower drag which accompanies laminar flow at 
high Reynolds numbers. 

Recently attention has been focused on the possible use of porous 
media as an enclosure for constant temperature environmental chambers. 

The problem of maintaining an enclosed area at a temperature within, such 
that it is suitable for human existance, when such an enclosure is sub- 
jected to extreme outside temperatures, has been encountered in the air- 
craft industry. Pilots who will experience atmospheric reentry frictional 
heating in the proposed new aircraft, must be provided with a cockpit which 
will remain at "shirt sleeve temperature". Hie use of porous walls has been 
proposed for the cockpit lining, with facilities for inducing air flow out- 
ward through the walls, thus cooling an otherwise radiant heat emitter. 

This system is believed to be an excellent solution to this problem. One 
primary advantage being that of great weight saving in air conditioning 
equipment, which would otherwise be required. 

These applications suggest the desireability of an increased program 
of study of the convective heat transfer behavior in porous solids. 

Although some work has been started in this country and more in Great 
Britian, there still remains a paucity of data useful for design at the 
present time. The principal reason for this lack, is the difficulty in 
accurate experimentation. The possibilities of using p>orous media for 
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transpiration cooling at low flow rates, provided the motive for the 
present investigation. 

i 
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2. Objectives* 

a. Develop a testing technique to determine the basic convective 
heat transfer data for the flow of air through wire screen matrices, at 
low Reynolds numbers, 

b. Obtain the heat transfer data for screen matrices over the 
Reynolds number range 0.1 to 10, based on hydraulic diameter. 

c. Summarize previous work in the field of porous media and to 
compare the previous results vdth the results of the present study. 

The work done in developing the analytical solution and method used 
in the present work is fully described by Locke in Reference (1)*. 



^Numbers in parentheses refer to Bibliography on page of this paper 
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3. Summary of Contents. 

A review of previous investigation procedures is presented first, 
followed by a discussion of the reasons this investigator chose the 
particular method used. The method of presenting the heat ''transfer data 
is then discussed. Presentation of the data will be made graphically as 
well as tabular ly. Graphic representation is plotted as Reynolds number, 
NRe’, as later defined, versus Ngt^Pr^^* The range of Reynolds numbers 
covered is from 0.1 to 10. 

A complete description of all the matrices tested is discussed. De- 
tails of each of the six types of woven wire screens are given, along with 
details of the matrices assembled from these screens. A complete descrip- 
tion of the test apparatus and test program is then presented. 

Heat transfer data for wire screen matrices, obtained by Coppage (2), 
are compared to the present results. Based on the results of this investi- 
gation recommendations for future work are made. 
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4. Review of Previous Work 



The theoretical analysis of the transient heat transfer behavior of 
a porous medium was first presented by Schumann (4) in 1929. The porous 
solid is considered to be homogeneous and initially at a uniform tempera- 
ture. A fluid at the same temperature as the solid, is flowing through the 
medium. At a certain instant the temperature of the entering fluid is 
assumed to change to a higher or lower temperature. The problem is then 
to find the temperature of the fluid and the solid as functions of time 
and position in the matrix. 

The analytical solution for the convective heat transfer coefficient 
of porous solids, requires that the following idealizations, initial condi- 
tions and boundary conditions be satisfied. 

Idealizations : 

1. The fluid specific heat and viscosity are constant 

2. The flow is constant 

3. The porous solid is homogenous 

4. The thermal conductivity of both the fluid and solid is 
infinite perpendicular to the flow and zero parallel to 
the flow. 

Initial and Boundary Conditions ; 

1. Initially the core is at a uniform temperature 

2. At time 0 = 0, the temperature of the incoming fluid changes 

instantaneously. ^ 

3. No heat passes the core boundaries. 

An energy balance on a fluid element within the porous medium yields 




(1) 
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The rate equation expressing the convective heat transfer between 
the porous medium and fluid is 



dft - X (± s - t f ) A oU (2) 

where X = total perimeter of the flow tube at a given 
cross section. 

An energy balance on the porous medium yields 

da, = - As f>i C s II s - d X (3) 

Introducing the variables '■£' and z for the variables 0 and x 
respectively, where 

- generalized time variable 

z = generalized position variable 

(JA)JL 

= NTU, when X = L 

and observing that the thermal capacity of the porous medium is much 
greater than the thermal capacity of the fluid within the porous 
medium, that is, W s c s )^WfCf, equations (1) and (2) combine to yield 



a** _ t - 1, 
s r 

and equations (2) and (3) combine to yield 

at e 



ar 



3 = t f - t £ 



The boundary and initial conditions are: 

tf(o,lT) = t fl 

tf (z,0) = tf + (tf^ - tf) e 



-z 



(4) 

(5) 

(6) 
(7) 
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The solutions to equations (4) and ( 5 ) together with the boundary 
conditions are found to be 
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( 8 ) 

(9) 



(10) 
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Schumann’s solution was first used as the basis for a transient testing 
technique by Furnas (5) in 1932 when he passed air through granular beds. 
The fluid temperature was measured at the exit of the matrix where z = NTU. 
Furnas compared his experimental data with Schumann’s theoretical constant 
z curves. The z value of the theoretical curve which fitted the experi- 
mental data was considered as the NTU of the matrix, and thereby the 
convection conductance h was evaluated. His work was limited in the 
Reynolds number range from 300 to 2000, based on the hydraulic diameter. 

Saunders and Ford (6) used essentially the same method in 1940 to 
obtain the heat transfer coefficient, h when air was passed through beds 
of spherical balls with porosity = 0.38. Tests were made at Reynolds 
numbers from 145 to 290, based on hydraulic diameter. Their results indi- 
cated a constant Stanton number, and with the assumption of Np r = 0.74 
for air, the following equation may be obtained from their work: 



N St N Pr 



2/3 _ 



0.05 
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Lof and Hawley (7) used the same transient technique in 1948 to ex- 
amine four sizes of gravel in beds over the Reynolds number range 5 to 500, 
based on hydraulic diameter. Their results indicated a large discrepancy 
existing between available data and even in the same group. The reason for 
this discrepancy, according to Coppage (2), may be due to the heat losses 
to the matrix boundaries and to other thermal capacitances in the experi- 
mental system. 

Romie (3) used the same transient technique in passing air through a 
screen matrix in 1948. He tested two wire screen matrices in which the 
screens were separated 1/32 inch, and covered the Reynolds number range 
from 20 to £00, based on hydraulic diameter and the mass velocity, based 
on free flow area. 

In evaluating the heat transfer coefficients, Romie plotted his data 
in the form (tf^ - tf^) / (tf^ - t^) versus 1 £'/N1U and measured the area 
under the curves up to a given value of ‘C/NTU. From Schumann* s theoreti- 
cal curves of (tf^ - tf) / (tf^ - ti) versus Y/z the area under the curves 
was obtained as a function of NTH, and the experimental NTU could thus be 
determined. 

Locke (1) in 1950 developed an improved method of evaluation of the 
experimental NTU from the theoretical curves to eleminate the possible 
error introduced by a displacement of the experimentally determined heat- 
ing curve either horizontally or vertically. In his method only the maxi- 
mum slopes of the theoretical and experimental curves need be compared. 

He derived an expression for the slope of the generalized heating curve, 

(tf - t^) / (t^ - t^) versus Y/z, finding that 




( 12 ) 
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From this it was 



for the position x=L, at which z = NTU and tf = tf . 

2 

found that the maximum slope of the generalized heating curve is a unique 
function of NTU, Therefore with this unique function available, only tho 
maximum slope of the experimental cooling curve need be determined. From 
the plot of maximum slope versus NTU, the NTU for the experimental curve is 
obtaiped and h is easily obtained from the relationship NTU=hA/wpCf. 

Saunders and Smoleniec (8) tested three screen matrices in 1951 in the 
Reynolds number range 20 to 1000, based on hydraulic diameter. They passed 
air through the matrix, and the convective heat transfer coefficient was 
evaluated by comparing the maximum slopes of the theoritical and experi- 
mental generalized curves. They presented the heat transfer data of screen 
matrices in a plot of Njj u versus Np e . When the curve was compared to the 
results of Njj u versus Nj^ e for flow normal to a single cylinder, the per- 
formance was quite similar. 

Coppage (2) tested six screen matrices in 1952 by using the same tran- 
sient technique developed by Locke (l). He covered the Reynolds number 
range 10 to 1000. Air was used as the working fluid. The porosities of 
the screen matrices tested were from 0.602 to 0.832, which covers almost 
the whole range of porosities of all square mesh, close packed screen 
matrices. In addition, he also tested a matrix of packed spheres. The 
diameter of the spheres was 0.0818 in., and the porosity of the sphere 
matrix was 0.390. 

This investigation has extended Coppage' s data to the Reynolds number 

2/3 

range 0.1 to 10. The data obtained is plotted as Nc t 'Np r versus Npj e ' 

and is compared to Coppage' s data at the overlap region of Np e ' approxi- 

2/3 

mately equal to 10. Coppage' s data is replotted as Ng^'Np^ ' vice 
Ngt'Npr -f° r comparison purposes. 
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5. Advantages and Disadvantages of Transient Method. 

The transient method as used by Furnas (5) and Saunders and Ford (6) 

is based on the solution of Schumann's (4) transient problem as outlined 

above. A porous solid initially at a uniform temperature is heated or 

tf 2 “ t>i 

cooled by a fluid passing through it, and — is plotted as a 



function of 
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this curve can be compared to the family of theoretical curves of (t^ - t^) / 
(tf^ - ti) versus TT/z with z as a parameter. Hie analytical curve which 
matches the experimental one will then give the value of N1U from which h 
can be computed. Hiere are three advantages to this method, (1) the fluid 
temperature only need be measured, (2) a core would be simple to construct 
and (3) any type of porous solid can be tested. A disadvantage of this 
comparison method is the difficulty encountered in the calculation of the 
theoretical curves. Calculation of these curves involves the evaluation 
of an infinite series for each point and the number of terms which muBt 
be taken to obtain a good approximation increases as NTU increases. 

As indicated above, Romie (3) obtained the area under his curves and 
compared them to the areas under the analytical solution curves. From the 
plot of the area as a function of NTU, the experimental NTU could be deter- 
mined and thus h could be obtained. As pointed out by Locke (l) both 
these methods have the disadvantage that a displacement of the experimen- 
tally determined heating curve either horizontally or vertically intro- 
duces errors in the results. 

In the method as developed by Locke (l) only the maximum slope of 
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the experimental plot of ( tf 2 - tj.) / (t fl - t±) versus T7NTU need be 
determined. Since the maximum slope of the generalized heating curve is 
a unique function of NTU, the experimental NTU value can easily be obtain- 
ed from the curve. This function from equation (12) is reproduced in 
Figure 8. A further advantage of the maximum slope technique is that the 
experimental range which can be conveniently covered is widened. In this 
method, actually, a premium is placed on a large NTU and IT in the analysis 
of the experimental data. Figure 9 is a plot of {% error in h) / (% error 
in max. slope) versus NTU. This plot shows the obvious advantage of a 
large NTU and suggests the avoidance of NTU in the range 1 thru 4 where 
the inflection point of equation (12) occurs. This did not pose a signifi- 
cant disadvantage in this investigation since porous media tend to have NTU 
greater than 5 because of large heat transfer area. In addition the flow 
rates used were sufficiently low in all cases to insure values higher than 
NTU =5. The chief disadvantage in this method is that any error in the 
experimentally determined maximum slope, is magnified in the corresponding 
value of NTU. However, a comparable disadvantage is also inherent in the 
other methods also. 

Thus for the usual porous media, with it's high NTU, the single blow 
transient method with the maximum slope technique is considered the most 
desirable and is used in the determination of data in this report. 

The method requires more care than the steady state technique to 
approach the idealizations and initial and boundary conditions of the 
theoretical analysis, which are restated below. 

Idealizations 

1. The fluid specific heat and viscosity are constant 

2. The flow is steady 

3. The porous solid is homogenous 
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4. The thermal conductivity of both fluid and solid is zero 
in the direction of flow and is infinite within the fluid 
and solid in the direction normal to the flow. Effective- 
ly all the convective resistance is lumped at the fluid- 
solid interface. 

Since the fluid used is air, the first idealization can be approached if 
temperature differences are kept small. For instance, if the maximum air 
temperature variation is from 80°F to 100°F, the maximum deviation from 
the arithmetic average will be plus or minus 1 . 5 % in viscosity and negli- 
gible in specific heat. 

The second and third idealizations can be approached by careful design 
and selection of equipment, and by discreet selection of the porous solid 
to be tested. The fourth idealization, however, can be approached only by 
correct core design. As pointed out by Locke (l) and Coppage (2), this 
poses no problem since heat transfer in the flow direction is ordinarily 
very small in the fluid and may be made negligible in the solid by inter- 
rupting the solid material. In the direction normal to the flow, the con- 
duction paths in both fluid and solid are usually small so that the effect 
is that of low thermal resistance within the solid relative to the resis- 
tance between them. 

Initial and Boundary Conditions 

1. Initially the core is at a uniform temperature. 

2. At time 9=0, the temperature of the fluid flowing pi st the 
screens changes instantaneously. 

3. No heat passes the core boundaries. 

The first condition is achieved satisfactorily if suitable precautions 
are taken to insure uniform heating of the screens at temperature other 
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than ambient. The second condition requires some method of rapidly intro- 
ducing the heated screens into the flow stream of constant temperature 
ambient air. The problem of heat loss to the walls of the approach ducting, 
as encountered by Coppage (2), is avoided in this particular technique, in 
that ambient air is flowing thru the test duct and thus remains constant up 
until the time it encounters the heated matrix. The third conditiai can 
only be minimized by careful design of the matrix holder. Since the screens 
are heated while contained in the matrix holding device, which has a certain 

thermal capacitance, an unavoidable transfer of heat from the holder to the 

$ 

screens will result as they are cooled. However this effect is minimized 
in the design of the matrix holder. 

Additional consideration of the foregoing requirements was made by 
Locke (1) and the design criteria established. 

Briefly, the experimental method used in obtaining the heat transfer 
data is as follows: the matrix, heated initially to a uniform temperature, 

is cooled by the fluid which enters at a constant, lower temperature. The 
temperature-time history of the fluid leaving the matrix is recorded dur- 
ing the cooling process. The maximum slope of the "cooling curve" is com- 
pared to that of the theoretical cooling curve. From a plot of the maxi- 
mum slope of the theoretical cooling curve versus NTU, the NTU correspond- 
ing to the maximum slope of the experimental cooling curve can be obtained. 
Thus, in tills manner the heat transfer unit conductance can be obtained. 
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6. Method of Presentation of Results. 

Through the dimensional analysis of the problem of heat transfer to 
a porous medium as developed by Coppage (2), it is seen that a definition 
for flow diameter, other than that normally used for flow thru pipes, will 
be necessary. In comparing two geometrically similar media it is clear 

i 

that only one diameter is needed to define the relative sizes of the flow 
passages. From the description of the porous media used in this investi- 
gation and previous investigations, it is seen that all such media have 
flow passages which are at least qualitatively similar, in that they all 
have a random, three-dimensional character and involve many contractions 
and expansions. Thus in comparing various porous media a single dimension 
can provide a very good measure of the relative sizes of the flow passages. 
One very logical choice of a dimension is based on the familiar concept of 
a hydraulic radius. The hydraulic radius of a prismatical flow passage is 
defined as 

- _ flow passage cross-sectional area 

h flow passage wetted perimeter"” 

If numerator and denominator are each multiplied by the length of the flow 
passage, then 

r _ flow passage volume 

h flow passage surface area 

This latter form now provides a possible basis for defining a suitable 
average hydraulic radius in a porous medium. Therefore, this definition 
will now be used for screen matrices even though the flow passages are not 
cylindrical. The relationship is valid since it is merely a defined quan- 
tity. 

In order to get the average r^ of all passages in the matrix, one can 
extend the above definition to 

— Total flow passage volume 

r h Total flow passage surface area 
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If both numerator and denominator are divided by the total volume of the 
matrix, the result is 

_ Total void volume/total matrix volume 

" Total void surface area/total matrix volume 



a 

" 0 

where 

ri - Total volume of voids _ pnpa^TTV 
Total volume of matrix 

/3 - Total surface area of voids 

^ Total volume of matrix 



The average flow passage diameter is the conventional hydraulic diameter 
D h = 4r h 

Nov; in order to establish an average mas;; velocity, an average free 
flow cross-sectional area, A^, may be defined such that 

A^L = Volume of flow passages in matrix 
= Volume of voids in matrix 

If both sides of the equation are divided by the total matrix volume we get 
Af = a At 



A^ ~ Total frontal area of the matrix 
Therefore, the average mass velocity 

G i= ^ s Jk G_ 

Af ocA t "a 

Likewise for the average flow velocity 

V 1 . W f w f V 

/®f A f C(p f A f ~ 

The above definitions are based upon the total matrix dimensions, and 
therefore some modifications must be made in order to adapt these defini- 
tions to a single screen of the matrix. To accomplish this, the following 
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definitions will be made 



& 1 = Porosity of a single screen 

ft « = Surface area / unit volume for a single screen 
Therefore 

r * CC 1 hydraulic radius of a single screen 
h = 0 ' 

A f » = cf'At 

G i' = w f / Af ' = G /cx' 

Vi* = w f / ^o f A f = V/^c ' 

Based upon these new definitions then, the dimensionless groupings for 
single screen matrices are: 

(h/ 0 1 'e f ) = 0 (Djj'Gi* j^f \ c f y^fAf iCC') s 0 (4r h ,G i'^fJ c f /<fAf ioC*) 
N S t’ = 0 (%e* > N Pr iOC') 

The heat transfer data obtained for the wire screen matrices are shown 
graphically in Figures 1 thru 6 and are plotted as N St ,N Pr 2 ^ 3 versus N Re* 
far various porosities. Each plot shows the results of each of the three 
screen matrices assembled from one type of screen. 



* 
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7. Description of Matrices. 

The screen matrix is packed in such a way that the individual screens 
are arranged in parallel and in close contact with each other. The direc- 
tion of the woven lattices of the screens are situated at random. The 
flow approaching and leaving the matrix is in the direction perpendicular 
to the screens. It is desirable that the test matrices be capable of ex- 
act geometric description, as regards surface area, volume of solids, 
diameters, etc., in order that the results may be put into forms which 
have general applicability. Consequently, matrices were chosen for which 
the geometry is completely controllable and yet which permit random three- 
dimensional flow which is characteristic of porous media. The reason for 
this particular choice is that wire screens permit a fairly wide varia- 
tion of geometric parameters while retaining essentially the same type of 
flow passage geometry. 

Test matrices were constructed from six different sizes of stainless 
steel wire screen. The table below summarizes the geometry of the test 
materials; the details of all matrices examined are given in Table I. 



Screen 


Nominal Mesh 
per inch 


Nominal Wire 
Dia., In. 


Screen ( 

Porosity * ct 


60 x 60 
-.011 


60 x 60 


.011 


.602 


60 x 60 
-.0075 


60 x 60 


.0075 


.675 


24, x 24 
-.014 


24 x 24 


.014 


.725 


16 x 16 
-.018 


16 x 16 


.018 


.766 


10 x 10 
-.025 


10 x 10 


.025 


.817 


5x5 


5x5 


.040 


.832 



-.040 
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The screen matrices vrere all composed of circular sections of stain- 
less steel wire screen having a diameter of six inches. For all the 
matrices tested each set consisted of three screens. Coppage's (2) data 
indicated that the number of screens had no affect on the convective heat 
transfer coefficient h. The prime consideration was to establish flow con- 
ditions which represented a true porous media. With such low flow rates 
(3 to 40 fv/hr), it was felt that these conditions were adequately met. 

All matrices were assembled in, and held in position during testing by a 
special matrix holding device described in Part III and shown in Figure 7. 

An important assumption made in the analysis for the transient temper- 
ature technique used in this investigation, is that of adiabatic boundary 
conditions existing at the matrix boundaries. It is exceedingly difficult 
to provide adiabatic boundary conditions in a transient situation, since 
all insulating materials have some capacity for storing thermal energy. 

By using balsa wood for the matrix boundaries and for the entire test 
section as well, the heat transferred from the boundaries was held to a 

y 

minimum. Locke (1) and Coppage (2) both point out the idealizations 
necessary in developing the experimental method used. From their analysis 
of heat loss to the boundaries they showed that the heat transferred from 
the boundaries to the screens was generally less than five percent of the 
storage capacity of the matrix up to the time of maximum slope. 
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8. Experimental Apparatus. 

The apparatus assembled for this investigation is entirely different 
from that used by both Locke (l) and Coppage (2) in their projects. A flow 
diagram of the test system is shown in Figure 10. Basically the system 
consists of a balsa test section with a sliding drawer for holding the test 
matrix, an air metering device, an air ejector to provide flow thru the 
test core, a fan type blower, an air heater, a flow measuring system, and 
a temperature measuring system. 

Ambient air is drawn through a 6 in. throat diameter, spun aluminum 
nozzle, by an air ejector system of a steam condensing plant. The nozzle 
was designed according to ASME specifications, contained in Reference (9), 
to insure smooth flow. The air then passes immediately into a 6 in. I. D. 
solid balsa test section. The test section is pictured in Figure 11. It 
is fitted with a sliding drawer which holds the three screen matrix. In 
all runs the orientation of the test apparatus is such that the incoming 
air flows vertically downward thru the nozzle and into the test section. 

The air passes thru the horizontal screen matrix and then past an air tem- 
perature thermocouple matrix ljjs in. downstream from the screens. From 
there it flows through a piping arrangement into one of three Fischer- 
Porter flowrators. Flow rate is controlled by a 3/4 in. needle valve on 
the vacuum side of the flowrator. Fine control is afforded by a £ in. 
by-pass needle valve around the large needle valve. With a maximum fluid 
velocity thru the nozzle of about 1 ft/sec, uniformity of pressure over the 
matrix section is insured. 

Adjacent to the test core is an air heating duct used for heating the 
screen matrix. The duct is an insulated L-shaped unit of circular cross- 
section, with a fan and heater located in the horizontal leg of the L. 
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This affords radiant shielding of the screen matrix thermocouples from the 
hot element of the heater. The duct is baffled in the vertical leg of the 
section of the unit in order to insure uniform air temperature discharging 
from the heater. The heating element and fan are individually controlled 
by two, 115 volt 60 cycle A.C., "Variacs" with a range of from 0 to 120 
volts. Because of the low flow rates used in the heating duct and with the 
baffles installed, very little adjustment of the variacs was found necessary 
to maintain constant temperature air. 

Temperature Measuring System . All thermocouples used are 24 gage 
Copper-Constantin. The flowrator temperature is measured by a mercury glass 
thermometer inserted downstream of the flowrator s. The upstream ambient air 
temperature is measured by two thermocouples located on the circumference 
of the bell of the inlet nozzle. The downstream air temperature was mea- 
sured by a matrix of five thermocouples. The entire downstream matrix con- 
sisted of thirteen thermocouples held in a frame and located in equal areas 
of the flow cross-section. A picture of the thermocouple matrix is shown 
in Figure 12. It was found that there was considerable temperature varia- 
tion near the walls of the test section. In constant temperature tests run 
on this thermocouple array, the five central thermocouDles, located in equal 
areas, gave excellent agreement. The outer two circles of thermocouples 
showed considerable variation, with the outermost circle giving the great- 
est deviation. The five central thermocouples in series, thus, were used, 
and measured the air temperature at the center of the core cross section. 

The reference junctions were all placed in an ice bath at 32°F. The 
leads from the downstream air temperature matrix were led to a Leeds & 
Northrop, "Speedomax", Model G, emf recorder with a variable range of from 
1 to 20 mv. The range selected for use on all runs was 0 - 2 rav. The 
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ambient air thermocouples were read with a Leeds & Northrup portable pre- 
cision potentiometer. Provision was also made to insure uniform tempera- 
ture of the screen matrix when heated in the air stream. Five thermo- 
couples were welded to the screen matrix with a pencil point welder. 

These thermocouples were also located in equal areas of the cross section. 
Reference for this system was also an ice bath at 32°F and were read with 

a Leeds and Northrup portable precision potentiometer. 

§ 

Sliding Drawer Matrix Holder . The experimental method requires that 
the test core be initially at a uniform temperature and that at the initial 
instant (©=0) air at a different temperature enter the upstream side. This 
sudden change in conditions is accomplished by sliding the balsa drawer 
matrix holder out of the balsa test section and over a heated air supply 
until the temperature of the screens is uniform and at the desired tempera- 
ture, The drawer is then pushed into the test section through which 
ambient air is flowing. 
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9, Experimental Procedure 

The conditions of the idealized problem upon which the analysis is based 
must be duplicated. Ambient air is drawn through the inlet flow nozzle and 
allowed to flow through the test matrix while the flow rate is adjusted. 

After constant flow is established the sliding drawer containing the screens 
is pulled out and into the heated air stream. When the screens reach a 
steady state temperature, as indicated by monitoring thermocouples spotted 
onto the screens, the drawer is quickly inserted into the air stream and 
the experiment is started. The ambient air temperature is recorded, the 
initial screen temperature is recorded, and the transient downstream air 
temperature is recorded on the Leeds & Northrup "Speedomax" recorder. 

After the downstream air temperature has returned to nearly ambient air tem- 
perature, the run is terminated. No detectible variation in flow rate could 
be noticed throughout the runs. The flow induced by imposing a vacuum on 
the downstream side of the test section, is seemingly, an excellent method 
to obtain low controlled flow rates. Considerable trouble was experienced 
with unsteady and erratic flow patterns when induced flow was attempted 
with the suction side of a centrifugal blower attached to the test section. 
The extremeiy steady flow rate was not obtained until the vacuum system of 
the air ejector steam condenser was used. 

The following data were recorded: 

a. Atmospheric pressure 

b. Atmospheric temperature 

c. Initial matrix temperature (5 Copper-Constantin thermocouples) 

d. Upstream ambient air temperature (2 Copper-Constantin thermo- 
couples) 

e. Downstream leaving air temperature as a function of time 

(5 Copper-Constantin thermocouples) 
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f. Flowrator temperature (Mercury thermometer) 
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10. Experimental Results . 

The geometric specifications of all the screen matrices examined are 
tabulated in Table I. All the test data for the matrices are evaluated 
on the basi3 of the geometry of an individual screen. In other words, the 
matrices are packed in such a way that the total length of the matrix is 
equal to the individual screen thickness multiplied by the number of 
screens. 

The experimental heat transfer performance is presented in Tables II 
to VII. The heat transfer data is also presented in graphical form by 
plotting versus N^ e ' for various porosities as shown in Fig- 

ures 1 through 6. 



26 



11. Discussion of Results 



As described by previous investigators, heat transfer behavior of por- 
ous media, appears to follow, essentially, a single characteristic. In 
looking at the plots in Figures 1 to 6, Coppage’s (2), data may quite sat- 
isfactorily be interpreted by straight lines drawn through the scatter. 

The general trend of the slope and position of the data is unmistakable. 
Slope and position of the lines show a definite increase with increasing 
porosity. 

Coppage (2) pointed out that some degree of curvature is suspected 
in his results. This statement is based on the fact that some of his data 
were taken in the transition range of Reynolds numbers. Friction factor 
data which Coppage took, showed a non-linear increase of friction factor 
with decreasing Reynolds number. Since the friction factor is a measure of 
form drag and skin friction drag, an increase in this combined coefficient, 
indicates a transition effect. Skin friction is affected very much less 
than form drag since the growth of even small standing vortices behind bluff 
bodies has a greater effect on form drag than on skin friction drag. There- 
fore, because of the analogy between heat transfer and skin friction drag, 
discussed in Reference (10), the heat transfer data would be expected to be 
more linear than the friction factor data but still exhibit some small cur- 
vature, Coppage did not predict the degree of curvature because of the 
lack of data, but rather used a straight line representation. Equipped 
now with additional data in the lower Reynolds number range, it can be 
seen that this curvature actually does exist. In the plots in Figures 1 to 
6, a proposed curvature is indicated which incorporates previous data with 
present data. 

In the extremely low Reynolds number range, it was noticed that a very 
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definite increase in the value of h occured. The explanation afforded 
this result is that at some low Reynolds number, the value dependent upon 
screen porosity, the effect of natural convective cooling affects the forc- 
ed convective cooling and causes unpredictable effects. The effects of 
this combined natural and forced convective cooling cannot be analyzed with 
the maximum slope technique used in this work. A technique for determining 
the natural convective cooling coefficient of porous media is certainly a 
problem which should be investigated. Since information is not available 
on this problem, no conclusions can be drawn from the results at these 
lowest Reynolds numbers. 

It should be noted that in comparing data obtained in the lower three 
screen porosities, w ith Coppage's (2) data, the curves are either smooth- 
ly connecting or overlapping with each other. Only on the highest three 
porosities was there considerable variation in results, and then only in 
the vertical position of the curve. The general slope of the comparable 
section of the curve is approximately the same. 

The data obtained will give only a slightly higher heat transfer 
coefficient h, than would result if the test section boundaries were tru- 
ly adiabatic. This is because of the heat which will flow into the 
screens from the boundaries during a run. In order to minimize this, the 
test core and matrix holder were constructed of balsa wood, a material 
with very low heat capacitance. Romie (3) and Coppage (2) point out that 
experiments support the statement that gain of heat from the walls, has 
little effect when measurements taken at the flow center-line are used 
and when the maximum slope technique is used in determining NlU. To fur- 
ther minimize the effects of non-adiabatic boundaries, both of these 
techniques are incorporated in the testing procedure used. 
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12. Comparison of Present Data with Previous Results . 

The heat transfer characteristics of square mesh screen matrices 
have been obtained over the Reynolds number range of 0.1 to 10. 

Heat transfer results for various matrices are correlated on the 
basis of the effect caused by porosity and Reynolds number. 

Previous data and present data agree quite well in their overlap 
region for the three lower porosity screens. There is a very definite 
curvature noted in the data as Reynolds number is reduced. There is no 
data available for comparison in the low Reynolds number range where this 
work. was done, however it does follow a curved pattern as predicted by 
Coppage (2). The amount of curvature is unknown, therefore only a propos- 
ed curvature is indicated. 

In the three highest porosity screens, vertical displacement of the 
curves occured. In the present data the value of h is somewhat higher 
than that shown by previous investigation. The possibility exists that 
in these higher porosity screens, a three screen matrix does not perform 
as a true porous media. It is possible, even with the very low flow rates, 
that there are insufficient obstructions to create the flow pattern re- 
quired of a porous media. If this is true, then a higher heat transfer 
coefficient will probably result for the lower flow rates. Coppage (2) 
stated that a lower coefficient for the first screen in a matrix should 
be expected, thus yielding a slightly lower overall coefficient. However 

4 

this appears to be true only at higher Reynolds numbers. At lower Rey- 
nolds numbers the opposite tendency seems to occur. Experiments with 
tube banks by Kays and Lo (11), show a definite tendency for h to be high- 
er for the first row of tubes than for the subsequent rows. Using this 
geometry for comparison then, it is entirely possible for the results to 
vary in the manner they do. 
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13. Conclusions. 

1. A testing technique has been developed for determining convec- 
tive heat transfer data for air flow through wire screen matrices at low 
Reynolds numbers. The method is also useful for testing any core with 
high effectiveness at low Reynolds numbers. 

2. Convective heat transfer data of square meshed, randomly packed 
screen matrices, composed of three screens, has been established for a 
range of porosities from 0,602 to 0.832 and for Reynolds numbers from 0,1 
to 10. 

3. Previous data for heat transfer to porous media are very limited 
in amount and non-existant over part of the Reynolds number range covered 
in this investigation. In the overlap region of data, excellent agreement 
was found in the three lowest porosity screens of the six sets of screens 
tested. However, further investigation is needed to substantiate the de- 
viation of the present data from previous data in the three higher porosity 
matrices. 

0 
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14. Recommendation for Future Work. 

a. No data is yet available for porous media on the heat transfer 
coefficient due to natural convection. A technique is needed for obtain- 
ing such data. 

b. Additional data is needed on the three matrices of highest poro- 
sity using more than three screens. This would verify whether the number 
of screens at the higher porosities have a significant effect on the heat 
transfer coefficient. 

c. In the screens of 5 x 5 - 0.041 in. mesh, it was not possible to 
obtain heat transfer coefficients for Np e * below 16. The erratic nature 
of the cooling curves was so severe that it was impossible to use them 

for data. Use of additional number of screens in the matrix should give 

* 

a much smoother response and enable extension of this matrix data to a 
much lower value of Np e ' • 

d. The analysis of the heat transfer coefficient of these screens 
by using the time temperature history of the screens, would afford an 
excellent check on the results of all the previous investigations. This 
could be accomplished by attaching a series of thermocouples to the screen 
wires by minute spot welding, and monitoring their time-temperature be- 
havior. Analysis of the data could be done by the thermal capacitance 
technique . 
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St'^Pr ' ver8u3 N^ e 1 for the 60 x 60 meah per inch — 0.011 inch wire diameter 

Figure 1 
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St'^Pr ^ ve r8U8 Npj e ' for the 60 x 60 mesh per inch —0.0075 inch wire diameter 
, Figure 2 
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versu3 K^ e 1 for the 16 x \C mesh per inch — C.C18 inch wire diameter 
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Photo shows the sliding drawer matrix holder. Ihe test 
screens are mounted in the drawer and held in place by a 
balsa frame. To the left of the drawer can be seen the test 
screens with thermocouples welded to them for measuring the 
initial screen temperature. 



Figure 7 
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Ambient Air Temp. 
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t ) — D esignates Temperature Measuring Points 




Photo shows the path traveled by air as it flows 
through and inlet nozzle and into the test core. 
Rectangular opening accommodates the sliding 
drawer which holds the screen matrix to be tested. 

Figure 11 
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Photo showing the thermocouple matrix consisting 
of thirteen thermocouples which could be read 
either individually or in series. Thermocouple 
leads were brought out of the edge of the drawer 
and into a selector switch. 

Figure 12 
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Photo shows the orientation of the test cell for 
all runs. The sliding drawer is resting over the 
warm air duct in the position required for heating 
the screens. The potentiometers were used to 
measure the ambient air and the initial screen 
temperature. The continuous pen recorder which 
was used to record the temperature-time history 
of the downstream air, is not shown. 

Figure 13 
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TABLE I 



DETAILS OF WIRE SCREENS USED IN MATRICES 



Screen 

Size 


60 x 60 
-0.011 


60 x 60 
-0.0075 


24 x 24 
-0.014 


16 x 16 
0.018 


10 x 10 
0.025 


5x5 

0.041 


Nominal 

Mesh/inch 


60 x 60 


60 x 60 


24 x 24 


16 x 16 


10 x 10 


5 x 5 


Nominal Wire 
Dia., In. 


0.011 


0.0075 


0.014 


0.018 


0.025 


0.041 


Actual Wire 
Dia., In. 


0.0106 


0.0075 


0.0136 


0.0177 


0.0245 


0.0405 


Frontal Area 
A t , ft 2 


0.1962 


0.1962 


0.1962 


0.1962 ^ 


0.1962 


0.1962 


Length of 
Matrix, ft 


0.00754 


0.00415 


0.00645 


0.00858 


0.0125 


0.0182 


Mass 

V , lhm 


0.123 


0.0653 


0.0873 


0.0962 


0.118 


0.150 


Porosity 
at ’ 


0.602 


0.675 


0.725 


0.766 


0.817 


0.832 


Heat Transfer 
Area, A', ft 2 


0.848 


0.582 


0.408 


0.343 


0.305 


0.254 


Heat Transfer 
area/unit vol. 
/ 3 ‘ , f t 2 /f t 3 


1815 


2090 


980 


623 


352 


196 


Hydraulic Radius 

r h ' , ft x 10 4 


3.32 


3.23 


7.40 


12.3 


23.2 


42.5 


Free Flow Area 

Af'sce'^t* ft 


0.118 


0.132 


0.142 


0.150 


0.160 


0.163 



Note: For all screens, material is 18/8 stainless steel, specific heat 

♦12 BIU/lbm°F} thermal conductivity, 10 BUJ/hr ft^ °F/ft} density, 
4 90 lhm/ft 3 . Parameters are computed for a single screen. 
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HEAT TRANSFER DATA 
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16 x 16 MESH PER INCH - 0.018 INCH WIRE DIAMETER SCREEN POROSITY - 0.766 
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0.550 2.80 99.88 65.95 36.25 4.64 0.400 



10 x 10 MESH PER INCH - 0.025 INCH WIRE DIAMETER SCREEN POROSITY - 0.816 
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